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ANALYSIS  OF  OSCILLATIONS  IN  THE  GYROTRON 
TRAVELLING  WAVE  AMPLIFIER 


I.  INTRODUCTION 

Current  experiments  performed  at  NRL  on  the  gyrotron  travelling  wave  amplifier  (gyro-TWA) 
have  yielded  record  power  levels  at  millimeter  wavelengths.' 1,2)  However,  the  presence  of  unwanted 
oscillations  has  thus  far  prevented  the  operations  of  the  amplifier  in  the  optimum  regime (2>.  Significant 
improvement  in  performance  is  anticipated  with  the  suppression  of  these  unwanted  oscillations.  The 
main  purpose  of  the  present  work  is  to  investigate  the  possible  causes  of  these  oscillations  and  thereby 
to  suggest  ways  for  their  suppression. 

Oscillations  are  generally  caused  by  either  an  external  or  an  internal  feedback  process.  The  exter¬ 
nal  feedback  is  provided  by  reflections  at  both  ends  of  the  interaction  region.  If  the  amplitude  of  the 
reflected  signal  reaches  a  certain  level  that  the  loop  gain  exceeds  unity,  the  wave  amplification  processes 
become  regenerative  and  oscillation  consequently  takes  place.  The  internal  feedback  is  a  result  of  the 
dispersiveness  of  the  unstable  medium.  Under  certain  conditions,  the  wave  may  grow  locally  without 
propagating  axially  out  of  the  system.  As  a  result,  large  amplitude  waves  can  simply  grow  from  noise 
level  perturbations.  Unlike  those  caused  by  end  reflections,  oscillations  produced  by  internal  feedback 
processes  received  comparatively  little  attention.  Their  considerations  have  special  importance  for  the 
gyro-TWA  in  that  the  optimum  operating  regime  of  a  gyro-TWA  happens  to  border  on  the  transition 
from  the  convective  instability  (amplification)  to  the  absolute  instability  (oscillation). (}>  Careful  study 
of  the  internal  feedback  processes  can  therefore  provide  a  guide  for  acheiving  optimum  operation  of  the 
gyro-TWA  without  the  excitation  of  absolute  instabilities. 

Experimentally,  there  is  convincing  evidence  that  both  types  of  oscillations  are  excited  in  the 

gyro-TWA  operated  at  NRL.  The  two  types  of  oscillations  have  different  characteristics.  Toward  the 
Manuscript  submitted  August  27,  1980. 
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end  of  this  paper,  we  shall  briefly  describe  the  experimental  results  in  the  light  of  the  present  theoreti¬ 
cal  study.  The  details  of  these  experiments  have  been  reported  elsewhere.  21 

We  shall  first  (Sec.  II]  present  a  quantitative  study  of  the  threshold  current  above  which  an  abso¬ 
lute  instability  would  be  excited  as  a  result  of  internal  feedback  processes.  The  influence  of  velocity 
spread  and  resistive  wall  loss  on  the  threshold  current  will  also  be  examined.  We  next  consider  oscilla¬ 
tions  caused  by  external  reflections  in  Section  III,  where  experimental  observations  of  both  types  of 
oscillations  are  reported.  The  major  results  are  summarized  in  the  abstract.  They  are  explained  in 
greater  detail  at  the  beginning  of  each  section. 

II.  INTERNAL  FEEDBACK  MECHANISMS-CONDITION  FOR 
THE  ONSET  OF  ABSOLUTE  INSTABILITY 

When  operated  near  the  'grazing  condition,”  the  mildly  relativistic  electron  beam  of  the  gyro- 
TWA  interacts  most  effectively  with  the  waveguide.14- 31  This  "grazing  condition'  is  achieved  when  the 
longitudinal  velocity  of  the  beam  is  equal  to  the  group  velocity  of  the  forward  wave  of  the  waveguide 
mode.  If  the  beam-wave  coupling  (determined  by  the  current  strength  and  magnetic  field  match,  etc.) 
is  sufficiently  weak,  the  beam  will  only  couple  to  the  forward  waveguide  mode.  The  gyro-TWA  would 
then  operate  according  to  what  is  expected  of  a  convective  instability. <3>  However,  if  the  coupling  is 
sufficiently  strong,  the  backward  wave  of  the  guide  mode  would  be  excited.  The  excitation  of  these 
backward  waves  might  then  provide  an  internal  feedback  mechanism,  resulting  in  an  absolute  instability 
which  would  manifest  itself  as  a  natural  oscillation  of  the  system.  (Under  these  conditions,  care  must 
be  exercised  to  intepret  the  gain  curve  of  the  "amplifier".'61) 

The  critical  current  which  marks  the  transition  to  absolute  instability  depends  sensitively  on  the 
magnetic  field.  It  is  obtained  by  using  the  stability  criterion  of  Briggs  and  Bers. 01  The  analysis  is  based 
on  the  assumption  that  the  system  is  infinitely  long.  The  characteristic  frequency  and  wavenumber  at 
the  transition  based  on  such  an  analysis  would,  however,  shed  light  on  the  required  length  (starting 
length)  of  the  system  if  it  is  to  avoid  these  natural  oscillations. 
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The  result  is  summarized  as  follows.  If  both  velocity  spread  and  wall  resistivity  are  absent,  and  if 
the  applied  magnetic  field  is  maintained  near  the  grazing  condition,  a  convective  instability  sets  in  at  a 
relatively  low  current.  This  convective  instability  becomes  an  absolute  one  if  the  current  (beam 
strength)  exceeds  a  certain  critical  value.  This  critical  value  is  given  analytically  [cf.  Eq.  (10)]  in  terms 
of  the  waveguide  dimension,  electron  parallel  and  perpendicular  velocities,  and  the  applied  magnetic 
field.  At  the  onset  of  the  absolute  instability,  the  oscillation  frequency  and  the  axial  wavenumber  are 
also  given  analytically.  It  is  found  that  typically  the  oscillation  frequencies  are  almost  identical  (but 
slightly  less  than)  with  the  cut-off  frequency  of  the  waveguide.  We  found  that  velocity  spread  in  the 
beam  has  virtually  no  effect  on  the  threshold  current  for  the  transition  to  absolute  instability.  The 
presence  of  a  small  amount  of  distributed  wall  resistivity  would,  however,  raise  the  threshold  current 
significantly.  It  can  be  shown  that  such  a  distributed  wall  loss  does  not  significantly  reduce  the  gain  of 
the  amplifier.'7’  In  this  regard,  the  presence  of  wail  resistivity  offers  the  best  hope  to  render  the 
amplifier  stable  at  high  current  level. 

In  Section  Il-a,  we  present  the  details  of  the  analysis  for  the  case  of  a  cold  beam  situated  in  a  loss¬ 
less  waveguide.  The  method  developed  there  has  been  extended  to  the  studies  of  a  beam  with  finite 
temperature  in  Section  Il-b.  Finally,  in  Section  II-c,  we  present  the  results  for  a  lossy  waveguide. 

a.  Cold  Beam,  Lossless  Waveguide: 

The  typical  configuration  of  a  gyro -TWA  (Fig.  1)  consists  of  an  annular  electron  beam  propagat¬ 
ing  inside  a  waveguide  of  circular  cross  section  of  radius  rw.  The  electrons,  guided  by  a  uniform  mag¬ 
netic  field  ( B0  e.),  move  along  helical  trajectories.  Ideally,  all  electrons  have  the  same  perpendicular 
velocity  vl0,  and  the  same  parallel  velocity  (v.0),  with  their  guiding  centers  uniformly  distributed  on  a 
surface  of  constant  radius  r0  (Fig.  1).  The  dispersion  relation  for  the  interaction  between  this  cold 
beam  and  the  TEm„  waveguide  mode  (m,  n  being,  respectively,  the  azimuth'll  and  radial  mode 
numbers)  at  the  s-th  cyclotron  harmonic  frequency  has  been  derived  in  Ref.  (S).  It  reads 
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_  -  4 vc2  f  (ft*2  ~  k2c2)  ( km„r0 ,  k„„rL0) 

Vo'iKm*  l  («  “  *:VjO“  5(lc)2 

(***  k.V.0)  Qsm  ^mn^O’  ^mn^LO^  1 

<o-k.vto- Sflc  Y 


(1) 


where  at  is  the  frequency  of  the  mode,  kz  is  the  axial  wavenumber,  y0  -  (1  —  v2Jc2  —  v; ?(/c2)"’/2, 
fto  ”  vi(/<\  nf  “  «*o/yo  «c,  rLo  -  v1(/nc,  kmn  -  *„„//•„,  AT**  is  the  n-th  root  of  J'm(x)  -  0.  •/*(*)  is 
the  Bessel  function  of  order  m,  J'm(x)  -  dJm(x)/<b r,  v  =  Ne2/ me2  is  a  dimensionless  beam  density 
parameter,  ZV  is  the  number  of  electrons  per  unit  axial  length,  and  the  functions  Kmn>  Q,m  are 
defined  as  follows: 


Kmn  -  J2(Xmn)  [1  -  m2/ X 2  „], 
Hm(x.y)  -  [y^UV’.O-)]2. 


and 


-  2 H,„(x,y)  +  y 


J}-  «  (x)J's(y)J"s(y)  +  „  _  ,(*> 


rt{y)r%.x{y)  -  |  Cy) 


Instability  occurs  at  wave  frequency  and  axial  wavenumber  such  that  a»2  -  k2c2  -  k2nc2  -0  and 
at  —  k.  v.q  —  jflf  ■■  0.  Hence  the  second  term  on  the  right  hand  member  of  Eq.  (1)  can  usually  be 
ignored  in  comparison  with  the  first  term.  We  can  further  approximate  at2  —  k2c2  by  k2„c2  in  the  right 
hand  side  of  (1).  Thus,  Eq.  (1)  is  simplified  to 

(at2  -  k2c2  -  k2nc2)  (at  -  k.v. 0-  silc)2  - - ~  (k2„c2)  0i2o  (2) 

yQrw  Kmn 

Henceforth,  we  shall  concentrate  on  the  stability  analysis  according  to  Eq.  (2). 


It  is  convenient  to  normalize  the  frequency  at  with  respect  to  the  cutoff  frequency  atc  =  kmnc  and 
k.  to  km„.  Then  Eq.  (2)  is  written  in  dimensionless  form  as 

(S2-  k2-  0(3-  it/3,,  -  b)2--t,  (3) 

where  S  =  at/atc\  k  =  kjkm„\  /9ti  =  v;0/c;  and  the  dimensionless  parameters  b  and  <  are  defined  by 


y^mXmKat2 


(4a) 

(4b) 
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Note  that  6  is  proportional  to  the  magnetic  field  and  c  is  a  measure  of  the  beam  strength.  Typically,  e 
is  a  very  small  number,  on  the  order  of  10-5  (or  even  less).  As  an  example*2’,  take  /n  —  0,  s  —  n  — 
1,  0iO  —  .4,  0iio  »  .266,  <uc/2 ir  —  34.3  GHz,  rw  —  .533  cm.  Then 

e  -  3.78  x  1(T6  /  (4c) 

where  /  is  the  beam  current  in  units  of  ampere. 

In  the  limit  of  zero  beam  strength  (e  »  0)  the  waveguide  modes 

£2  -  P  -  1  -  0,  (5a) 

and  the  beam  mode 

£  -  -6-0,  (5b) 

are  decoupled  [cf.  Fig.  (2)].  The  coupling  of  these  two  modes  is  strongest  when  the  dispersion  curve  of 

the  beam  mode  is  tangent  to  the  waveguide  mode  [cf.  Fig.  (2)].  This  occurs  when  the  magnetic  field  is 

adjusted  so  that  6  -  60,  where 

60=  (1-0,?)1/2.  (6) 

In  the  following,  we  shall  assume  that  0„  is  fixed,  and  that  6  is  in  the  vicinity  of  60.  We  shall  then 

study  the  properties  of  the  dispersion  relation  (3)  as  e ,  the  beam  strength,  is  varied.  This  way  of  fixing 

parameters  is  suitable  for  comparison  with  experiments. (1, 21 

We  state  here  the  principal  results.  The  details  of  derivations  of  these  results  are  given  in  the 
Appendices. 

(i)  If  6  >  60,  there  is  an  instability  whenever  c  >  0.  This  instability  can  either  be  convective  or 
absolute.  These  statements  are  proven  in  Appendix  A. 

(ii)  If  6  -  60  (i.e.,  exact  grazing  condition),  there  is  a  convective  instability  if  0  <  «  <  «<,.  This  ins¬ 
tability  becomes  an  absolute  one  when  «  >  «cg  where  eCK  is  defined  by 


c"4-  (27)l/4V^fc„. 

(7) 

“  («w  “  6o)  /  4  0n, 

(8) 

(60  +  6V2  0i?) 

(9) 

*  1  +  80,? 
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Furthermore,  at  the  onset  of  the  absolute  instability,  (i.e.,  e  —  «,*),  the  characteristic  wavenumber 
and  frequency  <LV  are  given  by  Eqs.  (8)  and  (9)  respectively.  Note  that  the  wavenumber  at  the  onset 
of  absolute  instability  is  purely  real.  This  general  result  is  demonstrated  in  Appendix  B.  The  proof  of 
the  rest  of  the  statements  in  this  paragraph  is  given  in  Appendix  C. 

(iii)  For  b  -  b0,  the  onset  of  the  absolute  instability  occurs  when  t  -  ee,  where  er  is  defined  by 


tj'4-  (27) 1/4  y/flii  k„ 

b  +  [8/8,?(l-62)  +  640,?] 1/2 


1  +  8/8,? 


(10) 

(ID 

(12) 


Again,  k,  and  5,  are  respectively  the  wavenumber  and  frequency  of  oscillation  at  the  onset  of  the  abso¬ 
lute  instability.  It  is  clear  that  Eqs.  (10)-(12)  reduce  to  (7)-(9)  when  b  -  b0,  i.e.,  when  the  gyro-TWA 
is  operated  exactly  at  the  grazing  condition. 


The  numerical  values  of  the  critical  beam  strength,  as  well  as  the  corresponding  frequency  of 
oscillations  and  wavenumber,  can  readily  be  computed  for  various  values  of  /3„  and  b  from  Eqs.  (7)- 
(12).  Listed  in  Table  1  are  the  values  generated  from  Eqs.  (7)-(9)  which  physically  correspond  to  the 
exact  grazing  condition.  Note  from  the  last  column  of  this  table  that  the  threshold  beam  current 
(measured  by  t„)  increases  with  /3„.  This  result  can  readily  be  threshold  beam  current  (measured  by 
«w)  increases  with  /3„.  This  result  can  readily  be  understood  by  noting  that  at  higher  axial  beam  velo¬ 
city,  disturbances  convect  away  at  a  higher  rate.  To  couple  the  forward  propagation  of  disturbances 
with  the  backward  waveguide  mode,  a  condition  for  the  excitation  of  absolute  instabilities ‘ 3 * ,  a  stronger 
current  is  required.  In  Table  2,  we  fix  /3„  at  .266.  We  then  calculate  the  threshold  current  at  various 
values  of  magnetic  field  strength.  We  note  that  if  the  magnetic  field  is  below  the  grazing  value,  a 
higher  beam  current  is  necessary  in  order  to  excite  the  absolute  instability,  because  the  beam  mode  and 
the  waveguide  mode  are  away  from  perfect  synchronism.  If  the  magnetic  field  is  above  the  grazing 
value,  the  unstable  spectrum  readily  extends  to  negative  values  of  k.  and  the  absolute  instability  is 
easily  excited.  In  the  event  that  the  beam  mode  (Fig.  2)  intersects  with  the  backward  waveguide  mode 
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the  absolute  instability  sets  in  as  long  as  the  beam  current  is  non  zero.  In  the  last  two  columns  of 
Table  2,  we  have  assumed  that  the  threshold  current  /  is  related  to  e  by  Eq.  (4c).  Note  further  that  the 
oscillation  frequency  of  these  oscillations  is  slightly  less  than  the  cutoff  frequency  of  the  waveguide. 

In  the  current  gyro-TWA  experiments  at  NRL11-2’,  /3„  —  .266,  <oc/2ir  —  34.3  GHz  and  kmn  -• 
7.18  cm~]  for  the  TE01  mode.  If  the  magnetic  field  is  maintained  at  the  grazing  condition  (b/b0  —  1) 
then  an  oscillation  at  a  frequency  of  34.26  GHz  is  predicted  to  occur  when  the  beam  current  is  raised 
beyond  .588  amps  as  a  result  of  absolute  instability  [cf.,  Row  (S)  in  Table  II].  Experimentally,  an  oscil¬ 
lations  at  34.27  GHz  was  indeed  observed18’  at  a  beam  current  of  1  amp  as  the  magnetic  field  is  raised 
to  the  vicinity  of  the  grazing  condition.  The  actual  threshold  current  observed  should  be  somewhat 
higher  than  the  predicted  value  as  the  latter  was  made  for  an  infinitely  long  waveguide  whereas  in  real¬ 
ity  a  finite  length  waveguide  was  used  and  coupling  out  the  ends  was  not  included  in  the  present  calcu¬ 
lation.  We  shall  discuss  in  greater  detail  the  experimental  results  in  Section  III. 

b.  Warm  Beam,  Lossless  Waveguide: 

We  have  studied  the  effects  of  the  velocity  spread  and  energy  spread  on  the  threshold  current 
above  which  an  absolute  instability  is  excited.  In  the  NRL  gyro-TWA  experiments,  the  energy  spread 
is  minimal,  of  order  1%  or  less.  The  spread  in  the  longitudinal  velocity  is  considerably  higher,  on  the 
order  of  10  per  cent.  We  found  that,  even  if  the  velocity  spread  is  as  high  as  20  percent,  it  has  a  negli¬ 
gible  influence  on  the  threshold  current.  In  other  words,  momentum  spread  cannot  stabilize  the  abso¬ 
lute  instability.  This  seemingly  surprising  result  can  be  understood  as  follows.  The  presence  of  an 
absolute  instability  requires  the  backward  waveguide  mode  to  be  coupled  to  the  forward  beam  mode13’. 
The  axial  wavenumber  k .  is  thus  constraint  to  a  value  in  the  vicinity  of  zero  at  the  transition.  [See  also 
column  3  of  Table  II].  Now,  the  velocity  spread  <  Av:>  modifies  the  dispersion  relation  by  an  amount 
on  the  order  of  (k.  <  Av.>2/ojc\  which  is  an  exceedingly  small  quantity  at  the  transition  to  the  absolute 
instability  by  virtue  of  small  value  of  k.,  leaving  the  threshold  current  calculated  in  the  previous  sub¬ 
section  essentially  unchanged. 
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The  above  conclusions  stem  from  our  analysis  of  the  dispersion  relation  which  includes  the  effect 
of  velocity  spread  and  energy  spread,  measured  respectively  by  <Av:>/v:0  and  <A> >/-y0.  Assuming 
that  these  quantities  are  small,  we  obtain  the  dispersion  relationship 


(S2  -  k2  -  1)1(5  -  ki 9„  -  b )2  -  \kp 


<Av.> 


_  <  A-y  > 

—  oj  - ‘ — 


This  equation  can  readily  be  compared  with  Gq.  (3),  where  the  symbols  are  defined.  Equation  (13)  can 
be  deduced  from  the  one  obtained  by  Uhm(”  who  applied  the  idealized  "water-bag"  electron  beam  dis¬ 


tribution  function  to  the  formulation  developed  earlier. <10) 


Applying  the  techniques  developed  in  the  Appendices,  we  found  that  a  velocity  spread  less  than 
twenty  percent  has  virtually  no  effect  on  the  transition  current,  for  reasons  explained  above.  An  energy 
spread  less  than  half  a  percent  actually  slightly  lowers  the  threshold  current. 

It  should  be  emphasized  that,  while  the  momentum  spread  is  not  effective  in  suppressing  the 
absolute  instability,  it  can  reduce  the  gain  and  bandwidth  of  the  amplifier  significantly'1*.  Once  the  fre¬ 
quency  of  the  input  signal  is  above  the  cutoff  frequency  of  the  waveguide,  k.  will  increase  significantly 
(in  comparison  with  a  value  close  to  zero).  The  velocity  spread  will  then  have  a  significant  effect  on 
the  gain  {cf.  Eq.  (13)1.  This  observation  is  corroborated  by  the  numerical  studies  on  the  gain  and 
bandwidth  of  the  amplifier  reported  in  Refs.  (1)  and  (7). 

c.  Cold  Beam,  Lossy  Waveguide: 

An  attempt  to  stabilize  the  gyro -TWA  is  to  coat  the  inner  wail  of  the  metallic  waveguide  with  a 
layer  of  material  of  finite  electrical  conductivity,  such  as  graphite.'8*  A  detail  study  of  the  effect  of  a 
resistive  wall  on  the  gain  of  the  gyro -TWA  is  beyond  the  scope  of  this  paper.  Here  we  concentrate 
only  on  the  dependence  of  threshold  current  on  wall  resistivity.  It  is  found  that  even  a  small  amount 
of  wall  loss  increases  the  threshold  current  substantially. 

Let  5  be  the  skin  depth  of  the  wall  material  of  the  waveguide.  Typically,  S/rw  «  1  where  rw  is 

the  wall  radius.  In  the  limit  of  small  8/r„„  the  dispersion  relation  reads'7* 
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SJ-  P-  [l“7“  (/-1)  |  ©-Wit 04) 
which  is  again  to  be  compared  with  Eq.  (3).  In  obtaining  Eq.  (14),  we  have  ignored  the  effect  of  velo¬ 
city  spread  as  it  can  hardly  modify  the  threshold  current.  In  passing,  we  remark  that  in  the  resistive 
wall  experiments  to  be  performed  at  NRL‘*\  8/r„  »  2.6  x  10“\  and  that  values  of  8 /rw  as  high  as  10~2 
are  reasonable  in  practice. 


The  analysis  of  Eq.  (14)  for  the  threshold  current  is  more  involved  and  will  not  be  detailed  here. 
Shown  in  Fig.  (3)  are  the  results.  To  obtain  the  data  in  this  figure,  we  fix  —  .266.  The  various 
curves  correspond  to  different  magnetic  fields.  They  illustrate  the  threshold  current  as  a  function  of 
8/rw  (wall  resistivity).  It  is  seen  that  a  small  amount  of  wall  resistivity  can  significantly  increase  the 
threshold  current.  That  is,  wall  resistivity  may  render  the  gyro-TWA  stable  against  the  absolute  insta¬ 
bility.  Shown  in  Table  III  are  the  data  for  b/  b„  —  .98  at  the  transition.  Again,  we  have  assumed  that 
the  beam  current  I  is  related  to  e  by  Eq.  (4c)  in  obtaining  these  results. 

The  increase  of  threshold  current  with  wall  resistivity  may  be  understood  qualitatively  as  follows. 
As  the  frequency  of  oscillation  at  transition  is  below  the  cutoff  frequency  of  the  waveguide,  the  back¬ 
ward  wave  of  the  waveguide  mode  attenuates  in  the  reverse  direction  of  beam  propagation.  This 
attenuation  rate  increases  as  wall  resistivity  (8 /rw)  increases.  For  an  absolute  instability  to  be  excited, 
this  attenuation  rate  must  match  the  amplification  rate  of  the  forward  beam  mode.  This  matching 
requirement  is  just  the  alternative  statement  that  a  saddle  point  (i.e.,  merging  of  roots  of  A;  at  a  single 
frequency)  must  be  present  foT  the  onset  of  absolute  instability. <M  A  higher  beam  current  is  thus 
needed  to  excite  the  absolute  instability  when  wall  resistivity  is  present.  Note  again  that  the  oscillation 
frequency  is  just  below  the  cutoff  frequency  of  the  waveguide  at  the  transition,  and  that  is  complex 
for  reasons  given  above. 

We  point  out  that  the  internal  feedback  mechanism  has  recently  been  theoretically  analyzed  by 
Etlicher  and  Buzzi.""  However,  the  main  stabilizing  influence  of  wall  resistivity  has  not  been  included 
in  their  calculations.  Their  work,  reported  in  a  conference  abstract'11*,  does  not  contain  sufficient  detail 
for  a  direct  comparison  with  our  work. 
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III.  EXTERNAL  FEEDBACK  MECHANISMS  AND  EXPERIMENTAL  OBSERVATIONS  OF 
OSCILLATIONS. 

Treating  the  length  of  the  amplifier  as  infinite  we  study  in  the  previous  section  the  convective  and 
absolute  nature  of  the  cyclotron  maser  instability.  When  an  absolute  instability  first  sets  in,  the  fre¬ 
quency  of  oscillations  is  predicted  to  be  just  slightly  less  than  the  cutoff  frequency  of  the  waveguide. 
The  axial  wavenumber  is  typically  1/30  the  transverse  wavenumber  of  the  waveguide.  Thus,  a  short 
interaction  length  (as  well  as  a  resistive  wall)  may  reduce  the  danger  of  exciting  absolute  instabilities. 

In  reality  the  gyro-TWA  system  is  finite  in  length  and  an  additional  mechanism  for  oscillation  is 
provided  by  partial  reflections  of  the  amplifying  waves  at  the  ends  of  the  system.  At  high  beam 
currents,  the  e-folding  length  of  spatial  amplification  is  very  short,  and  a  slight  mismatch  may  convert 
an  amplifier  into  an  oscillator.  Based  on  this  argument,  one  can  see  that  oscillations  due  to  end 
reflection  will  inevitably  set  in  at  a  high  current  level,  regardless  of  how  small  a  mismatch  there  is  at 
the  ends  of  the  device.  Oscillations  induced  by  end  reflections,  as  well  as  their  suppressions,  are  well- 
known  in  the  technology  of  travelling  wave  tubes.  For  the  present  case  of  the  gyro-TWA,  this  type  of 
oscillations  will  have  frequencies  higher  than  the  cut-off  frequency  of  the  waveguide  for  only  then 
would  the  gain  of  the  amplifier  be  significant.'1’  These  statements  appear  to  have  been  corroborated  by 
the  NRL  experiments  to  be  decribed  below.  Again,  a  resistive  wall  can  (1)  render  the  amplifier  stable 
against  this  type  of  oscillations,  and  at  the  same  time,  (2)  increase  the  bandwidth  of  the  amplifier, 
without  significantly  sacrificing  its  gain.'7’ 

Instabilities  have  been  routinely  observed'2,8’  during  experimental  operation  of  gyro-TWA 
operating  in  the  TE0l  circular-electric  mode  at  35.1  GHz  at  NRL.  The  physical  design  of  the  gyro-TWA 
uses  a  0.533  cm  radius,  for  a  cutoff  frequency  of  34.30  GHz,  copper  waveguide  21  cm  long.  The 
parameters  of  the  beam  were  designed  to  be  vn  -  .27  c  and  vx  -  .40  c.  Typical  linear  (small  signal) 
gains  obtained  were  24  db  at  a  3  ampere  beam  current  and  32  db  at  9  amps. 

Typical  oscillations  occur  at  several  discrete  frequencies  normally  just  above  cutoff  at  34.30  GHz, 
such  as  34.38  GHz  and  34.45  GHz.  The  power  output  of  these  oscillations  is  typically  on  the  order  of 
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the  maximum  saturated  output  during  stable  amplification.  These  oscillations  can  be  induced  from  a 
stably  operating  amplifier  by  increasing  the  beam  current  and/or  the  magnetic  field  hence  increasing  the 
gain.  These  oscillations  can  be  explained  as  follows. 

The  gyro-TWA  is  operated  near  the  cutoff  frequency  of  the  interaction  waveguide.  The  mechan¬ 
ism  actually  has  considerable  gain  almost  down  to  cutoff. 111  Even  though  the  input  and  output  couplers 
were  carefully  designed  for  low  reflections  throughout  the  useful  bandwidth  of  the  gyro-TWA,  match¬ 
ing  becomes  impossible  as  the  cutoff  frequency  is  approached.  Thus,  reflection  spikes  are  obtained  near 
cutoff.  Cold  test  reflection  measurements  indeed  reveal  reflections  at  the  problem  frequencies,  namely, 
at  34.38  GHz  and  34.45  GHz. 

As  the  gyro-TWA  beam  current  or  magnetic  held  is  increased  the  gain  near  cutoff  is  increased. 
The  reflections  constitute  a  feedback  and  when  the  net  loop  gain  exceeds  loop  losses  then  oscillation 
sets  in.  This  is  evidenced  by  the  oscillations  at  34.38  GHz  and  at  34.45  GHz  reported  above. 

Instability  at  just  below  empty  waveguide  cutoff  has  also  been  observed  and  is  believed  to  be  the 
absolute  instability  predicted  in  Sec.  II.  The  conditions  are  as  follows.  At  a  beam  current  of  1  amp  and 
cathode  voltage  of  -70  kV,  the  magnetic  held  was  scanned.  The  gain  measured  at  35.13  GHz  peaked  at 
12.7  dB  with  1230  watts  output  (unsaturated)  at  a  magnetic  held  of  12.81  kG.  As  the  magnetic  held 
was  increased  the  gain  (at  35.13  GHz)  dropped  rapidly  to  —0  db  at  13.11  kG.  Self-oscillation  in  the 
TEq{  mode  at  34.45  GHz  set  in  at  a  magnetic  Held  of  13.15  kG  at  a  power  output  level  of  8.71  kW. 
This  power  level  is  on  the  order  of  maximum  observed  saturated  amplifying  power  level.  An  observed 
(stable)  maximum  output  power  at  saturation  at  1  amp  beam  current  was  5.75  kW  with  4.3  db  gain.  As 
the  magnetic  held  was  raised  the  oscillations  at  34.45  GHz  continued.  At  13.37  kG  the  34.45  GHz 
oscillation  ceased  and  another  oscillation  at  34.27  GHz  was  observed.  This  oscillation,  slightly  below 
the  TEq\  cutoff,  had  an  output  power  of  200  watts.  At  13.41  kG  the  power  output  in  this  oscillation 
increased  to  660  watts. 
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The  experimental  data  reported  here  is  limited.  At  the  time  the  data  was  taken  n)  we  were 
unaware  of  potential  absolute  instabilities  or  their  characteristics.  Close  attention  will  be  made  to  abso¬ 
lute  instabilities  in  future  experiments.  However,  we  tentatively  interpret  the  above  data  as  follows. 

As  the  magnetic  held  is  increased  to  the  peak  gain  point  measured  at  3S.13  GHz  (corresponding 
to  a  grazing  condition),  the  gain  near  cutoff,  i.e.,  34.30  GHz  to  34.50  GHz,  is  too  low  to  sustain 
reflection  type  oscillations.  As  the  magnetic  field  is  increased  past  the  peak  gain  point  (at  35.13  GHz) 
the  gain  is  still  increasing  near  cutoff  until  sufficient  loop  gain  is  attained  for  oscillation  at  34.45  GHz. 
As  the  magnetic  field  is  increased  beyond  grazing,  the  threshold  current  for  the  absolute  instability 
rapidly  decreases.  When  the  absolute  instability  becomes  dominate,  then  oscillation  occurs  in  that 
mode.  The  calculations  of  the  threshold  indicate  that  the  absolute  instability  should  occur  even  at  graz¬ 
ing,  i.e.,  the  calculated  threshold  current  is  0.588  amp.  However,  the  calculation  was  made  for  an 
infinitely  long  waveguide  when  realistically  a  finite  length  waveguide  was  used  and  coupling  out  the 
ends  was  not  included.  Thus  we  would  expect  the  actual  threshold  to  be  somewhat  higher  than  indi¬ 
cated  in  Table  II.  The  lower  power  level  output  seems  reasonable  since  the  power  is  essentially  coupled 
out  the  end  rather  than  propagated  out. 

We  then  have  several  points  supporting  our  belief  in  the  observation  of  the  absolute  instability, 
(i)  The  oscillation  was  slightly  below  cutoff  in  TE^  mode  exactly  at  the  predicted  value;  (ii)  the  power 
output  level  was  substantially  lower  than  the  routinely  observed  reflective  oscillations  indicating  a 
different  type  of  instability;  and  (iii)  as  predicted,  the  absolute  instability  set  in  as  the  magnetic  field 
was  raised.  Although  the  reflection  type  instability  seems  to  be  the  most  dominant,  especially  at  higher 
current  levels,  this  absolute  instability  could  be  important  in  gyro-TWA  design  with  regard  to  methods 
of  overall  stabilization. 

An  attractive  method  to  suppress  both  types  of  oscillations  is  to  add  loss  to  the  circuit  to  decrease 
the  loop  gain.  This  also  raises  the  threshold  current  for  the  excitation  of  absolute  instabilities.  Further¬ 
more,  calculations'7’  show  that  the  forward  pin  in  a  distributed  loss  waveguide  circuit  is  only 
decreased  by  approximately  one-third  the  cold  tube  loss  in  dB  per  unit  length.  Thus  stability  would  be 
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obtained  with  only  a  slight  decrease  in  forward  gain.  In  addition,  wall  loss  leads  to  bandwidth  improve¬ 
ment. 
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APPENDIX  A 


Presence  of  an  Instability  Whenever  b  >  b0. 


In  this  Appendix,  we  prove  that  there  is  always  an  instability  (which  can  either  be  absolute  or 
convective)  whenever  the  magnetic  field  is  higher  than  that  at  the  grazing  condition.  It  suffices*3*  to 
show  that,  whenever  e  >  0  and  b  >  b0,  there  is  always  a  pair  of  complex  roots  of  a>  for  some  real 
values  of  k ,  say,  at  k  —  /9, Jb.  With  k  -  /3,,/A  Eq.  (3)  can  be  written  as 

O'2  -  lKy  -  a)2  -  -  e,,  (Al) 

where 

y  =  »b/(b2  +  pt)Ui,  (A2) 

a  9  (0J+  b2)1'2.  (A3) 

<|St  b*/(b 2  +  j8|f)2.  (A4) 

Note  that 


since 


a  >  l 


(A5) 


£,?  +  b:  -  1  +  ( b 2  -  bo)  >  I  if  b  >  b0. 


As  a  function  of  v,  /Cy)  5  (y2  -  DO'  -  a)2  is  sketched  in  Fig.  4.  As  easily  seen  from  this  figure,  the 
equation  /O')  “  -«i  admits  at  least  a  pair  of  complex  roots  of  y  whenever  <|  >  0.  From  (A2)  and 
(A4),  we  conclude  that  <o  has  at  least  a  pair  of  complex  roots  as  long  as  e  >  0. 
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Onset  of  Absolute  Instability  at  Real  Values  of  k. 

In  this  Appendix,  we  prove  that  when  an  absolute  instability  first  sets  in,  the  wavenumber  is 
always  real.  This  result  is  useful  because  it  considerably  restricts  our  search  in  wavenumber  space  in 
the  calculation  of  the  critical  beam  strength  that  marks  the  transition  to  absolute  instability. 


At  the  transition,  the  frequency  Z  -  u,  is  real.  For  this  value  of  Zs,  there  is  a  double  root13'  of 
k.  Below,  we  prove  that  this  double  root  is  real. 


I 

Since  Z  —  Sj  is  real,  and  Eq.  (3)  is  a  fourth  degree  polynomial  of  k  with  real  coefficients,  all  four 
roots  k  |,  k2,  ky,  k4  of  Eq.  (3)  must  either  be  real  or  appear  as  complex  conjugate  pairs.  We  shall  first 
assume  that  the  above-mentioned  double  root  has  a  nonvanishing  imaginary  part  and  then  show  that  a 
contradiction  results. 


Denote  the  real  and  imaginary  part  of  the  double  root  ku  k2  to  be  kv  and  k„.  If  ku  ^  0,  the 
other  two  roots  ky,  k4  must  also  be  a  double  root  since  the  roots  occur  in  conjugate  pairs.  In  this  case, 
Eq.  (3)  may  be  written  as 


(k  -  k\)(k  -  k2)(k  -  ky)(k  -  k4) 

—  (k  -  k„-  i  k„)2(k  -  k„  +  i  ksl )2  (Bl) 

-  k*  +  P  (—4 k„)  +  k2( 6 k2  +  2 k2) 

+  k  [-4 k„(k£  +  k,2)]  +  (kj  +  k2)2  -  0. 

On  the  other  hand,  equation  (3)  may  also  be  written  as 
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Since  (Bl)  and  (B2)  are  the  same  equation,  the  coefficients  must  be  the  same.  Thus, 


-  4  -  2(6  -  3j)/0n> 

\2 


6k,2.  +  Ik2 ■ 


b-  i 


fin 


-  (a*/  -  1), 


-  4 k„{k2  +  kj)  -  -2  (3, 2  -  1) 
Dividing  (B3)  by  (BS),  we  obtain 


6-3, 


fin 


—  (32  -  1)  -  kl  +■  k%. 

Using  (B6)  and  (B3),  we  express  the  right  hand  member  of  (B4)  as 


6  —  3, 


fin 


-  (o)}  -  1)  -  5 kl  +  k2,. 


(B3) 

(B4) 

(B5) 

(B6) 

(B7) 


Clearly,  (B7)  contradicts  (B4)  if  k,t  *  0.  We  conclude  then  that  the  assumption  that  ku  ^  0  is  false. 
Hence  at  the  transition  to  absolute  instability,  the  wavenumber  ks  must  be  real. 
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APPENDIX  C 


Detailed  Study  of  Absolute  and  Convective  Instability  at  the  Grazing  Condition. 

We  shall  focus  our  attention  to  the  case  where  b  -  b0.  The  dispersion  relation  (3)  reads 

Diot.k)  m  Ik2  -  (w2  -  1)]  k  -  ~~  -  -*7,  (Cl) 

011  011 

where,  in  this  Appendix,  we  have  dropped  the  bar  in  at  and  k  for  convenience.  Since  the  beam 
strength  e  enters  only  in  the  RHS  of  (Cl),  we  shall  examine  the  detailed  property  of  Diot.k).  In  the 
present  study  of  transition  to  absolute  instability,  we  need  only  to  consider  real  values  of  <0  and  k.  [Cf. 
Fig.  Si.  Since  the  onset  of  an  absolute  instability  always  occurs  at  a  multiple  root  of  real  k  (Cf.  Appen¬ 
dix  B),  we  need  only  to  focus  our  attention  at  those  values  of  k  in  Figs.  S  for  which  D(w,k)  is  station¬ 
ary.  Moreover,  since  «  ^  0,  only  non-negative  values  of  D(at,k)  are  considered. 

We  first  determine  the  trajectories  of  the  roots  in  the  complex  k  plane  as  the  imaginary  part  of  at 
changes  from  0  to  —  <».  (Cf.  Fig.  6a.]*  This  consideration  is  necessary  to  distinguish  an  absolute  from 
convective  instability.01  As  <0  —  -/  °°,  among  the  four  roots  of  k  in  Eq.  (Cl),  three  roots  would  tend 
to  — /«•  and  only  the  fourth  root  would  tend  to  +/00.  Note  that  the  fourth  root  (let  us  designate  it  to 
be  A 4)  corresponds  to  the  waveguide  mode  which  propagates  in  the  opposite  direction  of  the  beam  velo¬ 
city.  If  the  movement  of  the  roots  is  in  accordance  with  that  depicted  in  Fig.  (6b),  a  convective  insta¬ 
bility  sets  in  when  <  is  slightly  increased.  On  the  other  hand,  if  the  roots  move  according  to  Fig.  (6c), 
an  absolute  instability  would  set  in  when  c  is  raised  to  a  higher  level.  In  other  words,  an  absolute  insta¬ 
bility  sets  in  only  if  the  merging  of  roots  in  the  k  plane  involves  the  backward  wave  mode  of  the 
waveguide.  This  result  perhaps  is  anticipated  from  the  outset  because  the  internal  feedback,  which  is 
required  for  an  oscillation,  can  be  provided  for  only  by  the  wave  supported  by  the  waveguide  which 
propagates  in  the  reverse  direction  of  the  beam.  This  wave  can  be  excited  only  if  the  beam  strength  is 
so  high  that  all  modes  are  strongly  coupled. 

*We  assume  that  the  solutions  assume  a  dependence  of  exp  (mi  —  /Arc). 
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From  these  considerations,  we  then  see  that  the  merging  of  k  roots  at  k  -  0 tJb0  in  Fig.  (Sb)  must 


correspond  to  the  outset  of  a  convective  instability  since  the  merging  roots  arise  from  the  beam  mode 
and  the  forward  wave  mode  of  the  waveguide.  This  convective  instability  sets  in  as  long  as  e  >  0,  in 
agreement  with  the  conclusion  of  Appendix  A.  As  e  is  raised  slightly  above  zero,  we  encounter 
another  double  root  of  k  [cf.  kA  in  Fig.  (5c,d)].  This  double  root  of  k,  however,  does  not  correspond  to 
an  absolute  instability  since  the  merging  of  the  roots  does  rot  involve  the  backward  wave  of  the 
waveguide  mode  either. 


Examination  of  Fig.  5  shows  that  the  only  multiple  root  which  involves  both  the  backward  wave 
and  the  (forward)  beam  mode  is  kv  as  shown  in  Fig.  (SO.  Thus  the  onset  of  absolute  instability 
occurs  when  e  -  ecg,  where  tclt  is  given  by  {cf.  (Cl)l 


and  wv,kv  are  the  real  root  determined  from  the  conditions 


dD 

dk 

b1D 
8  2k 


-  0. 


-  0. 


A  little  algebra  yields  the  relevant  solution 


(C2) 


(C3) 

(C4) 


bo  +  6>/20i? 


w  1  +  80,? 

“  7T"  -  *o)- 

ii 


(C5) 


(C6) 


These  equations  are  just  in  Eqs.  (9),  (8)  of  the  main  text.  Substitution  of  (CS)  and  (C6)  into  (C2) 
gives  Eq.  (7). 


A  few  remarks  concerning  our  calculations  are  in  order.  It  can  readily  be  shown  from  Eq.  (9)  that 
is  always  less  than  unity.  This  fact  that  oscillation  may  occur  at  a  frequency  lower  than  the  cutoff 
frequency  of  the  waveguide  is  hardly  a  surprising  result.  It  is  well-known"2’  that  when  a  waveguide  is 
loaded  with  a  plasma,  modes  with  a  frequency  below  the  cutoff  frequency  of  the  waveguide  can  exist. 


LAU,  CHU.  BARNETT,  AND  GRANATSTEIN 


In  the  present  case,  the  waveguide  is  partially  loaded  with  the  beam,  which  itself  constitutes  an  active 
medium.  This  brings  out  another  point  which  at  first  sight  may  also  seem  a  bit  paradoxical.  At  the 
grazing  condition,  one  may  then  ask,  how  can  an  absolute  instability  ever  exist  as  the  beam  mode  and 
the  waveguide  mode  both  propagate  in  the  forward  direction?  The  answer  is  that  at  very  low  beam 
current,  both  modes  are  weakly  coupled  and  indeed  no  absolute  instability  exists  according  to  the  usual 
notion  of  coupling  of  modes.  However,  as  a  sufficiently  high  current  level,  both  modes  are  strongly 
coupled.  In  fact,  all  modes  are  coupled,  including  the  waveguide  mode  which  propagates  in  the  oppo¬ 
site  direction  of  the  beam.  The  excitation  of  the  latter  modes  provides  an  internal  feedback  and  is  the 
basic  reason  for  the  possible  excitation  of  an  absolute  instability.  [C f.  Fig.  (50  J. 
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LAU.  CHU,  BARNETT.  AND  GRANATSTEIN 


Fij,  3  —  Threshold  current  (1)  as  a  function  of  skin  depth  (wall  resistivity)  at  various  mainetic  field  streit|th  (b/b„). 

Here /3 1 1  —  0.266,  wr/2ir  -  34.3 GHz  and  r,  -  0.533  cm. 
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the  function  f(y)  -  O'1  -  l)(y  -  a)1 


LAU,  CHU,  BARNETT,  AND  GRANATSTEIN 


Fig.  6  —  Trajectories  of  the  roots  in  the  k  plane  as  w  changes  from  <•>,  to  u,  -  /«> .  In  (b),  a  convective  instability  exists. 

In  (c)  k,  corresponds  to  the  onset  of  an  absolute  instability. 
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